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We update the study of the total if/ and T production cross section in proton-proton collisions at RHIC ener- 
gies using the QCD-based Color-Singlet (CS) Model, including next-to-leading order partonic matrix elements. 
We also include charm-quark initiated processes which appear at leading order in a s , but which have so far been 
overlooked in such studies. Contrary to earlier claims, we show that the CS yield is consistent with measure- 
ments over a broad range of J/if/ rapidities. We also find that charm-quark initiated processes, including both 
intrinsic and sea-like charm components, typically contribute at least 20 % of the direct J/ \fr yield, improving 
the agreement with data both for the integrated cross section and its rapidity dependence. The key signature for 
such processes is the observation of a charm-quark jet opposite in azimuthal angle <p to the detected J/tf/. Our 
results have impact on the proper interpretation of heavy-quarkonium production in heavy-ion collisions and its 
use as a probe for the quark-gluon plasma. 
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The hadroproduction of J /if/ and Y is one of the key topics 
in phenomenological QCD. As opposed to lighter mesons, it 
is a priori straightforward to compute their production rates 
from gluon-induced subprocesses such as gg — > Qg (Fig. 1 
(a)), particularly since one can use nonrelativistic approxima- 
tions. However, there are many outstanding theoretical issues, 
including the role of color-octet (CO) states, the impact of 
next-to-leading order (NLO) - and even higher order - QCD 
corrections (Fig. 1 (c,d)), and the role of hard subprocesses 
such as gc — > J/ific (Fig. 1 (b)) which utilize the c-quark distri- 
bution in the proton. Other issues include the J/ ifi polarization 
puzzle, the factorization-breaking strong nuclear dependence 
in J/ifj hadroproduction at high Xf, and the uncertain effects of 
rescattering and energy loss mechanisms. All of these issues 
have impact on the proper interpretation of heavy-quarkonium 
production in heavy-ion collisions and its use as a probe for 
the quark-gluon plasma. For recent reviews, see [1]. 

It is widely accepted that a\ and a 5 s corrections to the 
CSM [2] are fundamental for understanding the pj spectrum 
of J/ifi and Y produced in high-energy hadron collisions [3- 
8]. However, if anomalously large contributions to the total 
cross section arise from at graphs, this would cast doubt on 
the convergence of the expansion in a s . It is thus important to 
check that LO and NLO predictions are close to each other and 
in agreement with experimental data. In this paper we carry 
out the first theoretical analysis in the CSM at NLO accuracy 
of the total J /iff, ifr(2S), and Y production in pp collisions at 
the BNL RHIC. We show that hard subprocesses based on 
CS QQ configurations alone are sufficient to account for the 
observed magnitude of the /^-integrated cross section. In par- 
ticular, the predictions at LO [2] and NLO [3, 4] accuracy are 
both compatible with measurements by the PHENIX collabo- 
ration at RHIC [9] within present errors. We shall also show 
that hard subprocesses involving the charm quark distribution 
of the colliding protons (Fig. 1 (b)) which constitute part of 
the LO (a^) rate, are responsible for a significant fraction of 
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FIG. 1: Representative diagrams contributing to 3 5 1 quarkonium (de- 
noted <2) hadroproduction via CS channels at orders al (a,b), 
(c,d,e,f). The quark and antiquark attached to the ellipsis are taken 
as on-shell and their relative velocity v is set to zero. 



the observed yield. Reactions such as gc — > J/if/c (thereafter 
referred to as eg fusion) also produce a charm jet opposite in 
azimuthal angle to the J ft//; furthermore, the rapidity depen- 
dence of this "away-side" correlation is strongly sensitive to 
the mechanism for the creation of the c-quark in the proton. 
An analysis of the invariant mass distribution of the J/tf/ + D 
pair may also shed light on possible contributions beyond the 
color singlet model, as described by the Color Transfer Mech- 
anism (CTM) [10, 11]. 

Subprocesses involving eg fusion with a charm quark from 
the proton have been considered in [12, 13] with the main fo- 
cus on the high pj spectrum At low pj, the typical scale of the 
production process is rather small, and thus one does not ex- 
pect higher-order QCD corrections such as gluon splitting into 
cc to give a significant contribution to the total cross section 
For example, the contribution to the total cross section from 
the process gg — > J/if/cc, appearing at a A s (Fig. 1 (e)) [14], 
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is at the level of 0.5 %. In contrast, in the case of intrinsic 
charm (IC) contributions, the c and c quarks are created from 
two soft gluons connecting to different valence quarks in the 
proton as in the BHPS model [15]; such contributions are rele- 
vant to charmonium production at all scales. The contribution 
from cy fusion was studied in photoproduction in [16]. 

We shall focus here on the "direct" hadroproduction of the 
J/if/, i//(2S), and Y(15) without the contribution arising from 
the decay of heavier states; this avoids the discussion of the 
production mechanisms of f-waves which are not well un- 
derstood. Although the total cross section for L — 1 states 
has been studied at NLO [17], an effective evaluation of the 
production cross section requires the introduction of an in- 
frared cut-off (as for their decay [18]) or CO contributions [19] 
which introduce new unknown non-perturbative parameters. 
Furthermore, the impact of the off-shellness of initial gluon 
on the^ c i yield may be significant [20, 21]. We have also re- 
stricted our analysis to the integrated-p^ distribution. Indeed, 
as noticed at the Tevatron energy [3, 4], the NLO pj distri- 
bution, contrary to the integrated one, can be negative at low 
pj. In addition, initial-state radiation [22] would also be ex- 
pected to significantly modify the spectrum at small pj and to 
increase (pi). 

In the case of Jj^i hadroproduction, the PHENIX data [9] 
includes the direct yield, but also a B feed-down fraction 
(4 + _ 3 2 % [23]), a if/(2S) feed-down (8.6+2.5% for \y\ < 0.35) and 
ax c feed-down estimated to be < 42% at 90% C.L. [23]. A re- 
cent analysis [25] from fixed-target measurements in pA sug- 
gests that it amounts to 25 + 5%, while the CDF measurement 
in pp at Fermilab gives 30 ±6% of the prompt yield for pj > 4 
GeV [26]. For our analysis, we will make the hypothesis that 
the^f c feed-down fraction is 30± 10% of the prompt yield inde- 
pendent of rapidity. Overall, we shall take F^ ct = 59 ± 10% 
and multiply the PHENIX results by this factor. The differ- 
ential J/ip production cross section vs y has been measured 
by PHENIX in the central (|y| < 0.35) as well as in the for- 
ward (1.2 < \y\ < 2.2) regions [9, 24]. The extrapolation to 
the direct yield using FjJ^ ct — 59 ± 10% is shown on Fig. 2 
(a). For the i[/(2S), only a negligible B feed-down competes 
with the direct mechanism. The preliminary measurement by 
PHENIX is shown on Fig. 2 (b). The T(15 +2S + 3S) cross 
section has been measured by STAR [27] and PHENIX [24] 
in the central region, and by PHENIX [28] in the forward re- 
gions. From the CDF analysis [29] at p T > 8 GeV, 50% of 
the T(15) are expected to be direct. Using the relative yields 
from [30], we expect 42 ± 10% of the T(15 +25 +3S)signal 
to be direct T(15). PHENIX and STAR data mutiplied by this 
fraction are displayed on Fig. 2 (c). 
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where P — pq + pg, p = (pq - pg)/2, s\,S2 are the 
heavy-quark spin and 5" / -\fN~c is the projector onto a CS 
state. In the non-relativistic limit, N(A\s\,S2) can be written 

as ^| — v(|, S2)yu(j, si) where is the polarisation vec- 
tor of the quarkonium. The sum over the spins yields to traces 
evaluated in a standard way. 

In our evaluation, we use the partonic matrix elements from 
Campbell, Maltoni and Tramontano [3] to compute the LO 
and NLO cross sections from gluon-gluon and light-quark 
gluon fusion. We guide the reader to [3] for details concerning 
the derivation of M(ab — > Q A (P) + j) at a s , the corresponding 
expressions at or 3 can be found in [33]. In the case of the eg fu- 
sion (at LO), we use the framework described in [34] based on 
the tree-level matrix element generator MADONIA [35]. For 
the parameters entering the cross section evaluation, we have 
taken \Rj^(0)\ 2 = L01 GeV 3 and |^^ (2S) (0)| 2 = 0.639 GeV 3 . 
We also take Br(//^ -> l + (T) = 0.0594 and Br(^(25) -> 
( + (-) = 0.0075. For the T(15), we will choose |/?(0)| 2 = 
7.6 GeV 3 , and Br(T -> ft) = 0.0218. The uncertainty 
bands for the resulting predictions are obtained from the com- 
bined variations of the heavy-quark mass within the ranges 
m c = 1.5+0.1 GeV and m b = 4.75 ± 0.25 GeV, 1 the fac- 
torization p F and the renormalization p R scales chosen 2 in 
the couples ((0.75,0.75); (1, 1);(1,2);(2, 1);(2,2)) x m T with 
m\ = 4m,Q + p T . Neglecting relativistic corrections, one has 
in the CSM, Mj/y = M^s) = 2m c and M r = 2m b . The 
parton distribution used was the LO set CTEQ6_L [38] for 
the LO gg fusion, the NLO set CTEQ6JVI for the gg + gq 
NLO one and, for the eg fusion, the LO set CTEQ6.5c [39] 
based on a recent global PDF fit including IC. We have em- 
ployed three choices for the charm distribution: (i) without 
IC [c(x,po) = (po =1.2 GeV)], (ii) with BHPS IC [15] 
((x) c+ c= J x[c(x) + c(x)]dx -2%) and (iii) with sea-like IC 
((x) c+ c = 2.4%) . While there does exist an intrinsic fr-quark 
content in the proton scaled by m 2 /;«^ relative to IC, its cor- 
responding contribution to T + b is additionally suppressed at 
RHIC energy by phase space due to the presence of an addi- 
tional fo-quark in the final state. 

We now describe our results. As shown in Fig. 2 (a) and (b), 
the yields at LO and NLO accuracy are consistent in size, and 
the uncertainty of the latter one (indicated by the two curves 
in both cases) is smaller than that of the LO. This provides 



In the CSM [2], the matrix element to create a 3 S i quarko- 
nium Q of momentum P and polarisation A accompanied by 
other partons, noted j, is the product of the amplitude to 
create the corresponding heavy-quark pair, a spin projector 
N(A\si, Sz) and R(0), the radial wave function at the origin in 
the configuration space, obtained from the leptonic width [32], 



It is common to see a wider range used for m c in NLO evaluations of open- 
charm cross sections, i.e. m c = 1.5±0.2 GeV (see [36]). In the case quarko- 
nium production within the CSM, such values so different from Mq/2 may 
require the inclusion of non-static corrections, which is beyond the scope 
of our analysis. See also our comment regarding the \j/(2S ) results. 
2 In principle, the renormalization scale ambiguity can be removed using the 
method described in [37]. 
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FIG. 2: (a) dcrff^'/dy x Br from the measurements by PHENIX [9, 
24] multiplied by our estimate of F^ ct compared to the CSM at LO 
(a 3 s ) by gg fusion only (thin-dashed lines), at NLO (up to aj) by 
gg and qg fusion only (thick-solid lines) and the sum "NLO + eg 
fusion" with the sea-like c(x) [39], denoted NLO + (light-blue band), 
(b) same as (a) for the iff(2S) with PHENIX data [24]. (c) same as (a) 
for the direct T with STAR [27] and PHENIX [24, 28] preliminary 
measurements for Y(15 +25 +35) multiplied by our estimate of 
F!j,' rccl (without NLO + , see text). The gaps between the two solid and 
the two dashed lines as well as the band reflect the variation of the 
cross section after a combined variation of the scales and the masses 
as indicated in the text. 



The yields at LO and NLO accuracy are compatible with the 
PHENIX data, in contrast to the conclusion of [40], in which 
feed-down from XcO an d X& at a 2 s was incorrectly assumed 
to be the dominant source of J/if/ production. This supports 
the good description of STAR results [42] for the J /if/ differ- 
ential cross section at mid pj predicted by the CSM at NLO 
including leading-/?/- a\ contributions (NNLO*)[5]. Note that 
a significantly larger CS yield points to a small impact from 
s -channel cut contributions [41]. 

Even though the NLO is close to the data, the additional eg 
contribution (even with a sea-like IC distribution) improves 
the agreement. However, phase-space effects are not properly 
taken into account in the case of iff(2S) production due to the 
restriction M^gsy = 2m c . The i//(2S) case is nevertheless en- 
couraging, since it does not involve the uncertainties arising 
from the extrapolation of the experimental data to the direct 
yield. We also give in Fig. 3 our prediction at tJs = 500 GeV 
for the direct J /if/ and Y yield for future comparison with the 
data taken this year. 
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FIG. 3: dcrfdy x Br for the direct yield of J /if/ and T as a function of 
y at i/S = 500 GeV for the same parameter ranges as Figs. I 



We note that the contribution from eg fusion (the results la- 
beled 3 NLO + were obtained with the sea-like IC from CTEQ 
6.5c) is significant for both J fx// and if/(2S) production and 
calls for a deeper analysis. First, It should be noted that 
NRQCD factorization breaking effects, such as those arising 
from the CTM [10, ll] may impact the low pj region. Such 
effects arise from infrared sensitive domains at NNLO when 
the 3 heavy quarks have comparable velocities. A careful 
study of the CTM is however beyond the scope of our analy- 
sis. Second, to precisely assess the impact of other choices for 
the charm distribution, c(x), we have evaluated the fraction of 
J I if/ produced in association with a single c-quark relative to 



3 Notation not to be confused with NLO* or NNLO* which denote real- 

some indication that we are in a proper perturbative regime. emission contributions as evaluated in [5]. 
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FIG. 4: Fraction of J /if/ produced in association with a single c-quark 
(via gc — > J/if/c) relative to the direct yield (NLO + ) as a function of 
y$ and for three models for c(x): without IC (No IC), sea-like and 
BHPS (see text). 



the direct yield as a function of and for the three models 
for c(x). These are displayed on Fig. 4 for which we have set 4 
m c = 1.4 GeV and varied pp and pa within the same values as 
for Figs. 2. This clearly confirms the impact of the eg contri- 
bution, which ranges from 10 % up to 45% of the direct yield 
in the case of sea-like c(x). 

Note also that at larger pj, we expect significant a s 
contributions from eg fusion, since they then exhibit a 
fragmentation-like topology (Fig. 1 (f)). This was studied by 
Qiao [13] for the Tevatron using a conventional c-quark dis- 
tribution, but this evaluation cannot be extended to small pj 
where it is infrared divergent. For the BHPS IC distribution, 
the pj distribution at large p T and RHIC energy will show 
an analogous enhancement as seen at large rapidity in Fig. 4. 
This may also impact the J /if/ yield in this region. In order 
to assess experimentally the importance of eg fusion, whether 
from the usual CSM or from CTM effects, the measurement 
of J /if/ in association with D meson would be illuminating, as 
has been noted in ref. [14] for J /if/ + cc. More accessible is 
the study of the azimuthal correlation of J /if/ + e in the central 
region by PHENIX and STAR and of J /if/ + p. in the forward 
region by PHENIX. The key signature for such subprocesses 
is the observation a lepton excess opposite in azimuthal angle 
<f> to the detected J /if/. 

One can also have, at large rapidity, (cc)g — > J /if/ contribu- 
tions to the total cross section [44, 45] from the coalescence 
of the charm pair and gluon; in this case the J/if/ acquires the 
momentum of both the c and c quarks from the projectile. In- 
trinsic charm Fock states such as \(cc)% c {uud\ c ) can explain 
J/i/j and double J/ifi production at high Xf > 0.6 observed in 
pA and nA collisions by the CERN NA3 experiment as well 



4 Indeed, for our prediction of the ratio to make sense, the colour singlet 
contribution has to be the dominant one, which can onlybe the case for a 
rather low charm quark mass such as m c = 1 .4 GeV. 



and it anomalous A dependence [31]. 

We now turn to T hadroproduction where the bg fusion pro- 
cesses are suppressed by phase-space and by the 1 /m? depen- 
dence of the b-quark content in the proton. Thus we have only 
computed the LO and NLO yield from gg and qg (see Fig. 2 
(c)). The predictions are not far from the extrapolation of pre- 
liminary data by PHENIX and STAR. In addition, the con- 
sistency between CDF data at the Tevatron at mid and large 
pj and the very first NNLO* CS analysis [5] also suggests 
that T production can be understood from perturbative QCD. 
We also emphasize here that the rapidity region accessible at 
RHIC allows for measurements of T production at high xp 
very close to 1 where the intrinsic bottom quark pair can sim- 
ply coalesce to form a T after a single scattering to change its 
color in {bb\ c + g — > T in analogy to the large Xf J /if/ pro- 
duction [45]. It does not require a third b-quark and is thus 
not suppressed by phase-space effects. 

We now briefly discuss the production of Jfi/i in pA colli- 
sions as CS states, likely the dominant mechanism at RHIC 
energy. In the central region, the cc pair hadronizes outside 
the nucleus. Although the energy loss of a colored object 
in cold nuclear matter is limited to be constant, rather than 
scaling with energy, by the Landau-Pomeranchuk-Migdal ef- 
fect [48], its magnitude per unit of length will be significantly 
larger for a CO than for a CS state. The recent observation 
by STAR [42] of the non-suppression of J/if/ in Cu-Cu col- 
lisions at increasing pj clearly supports the hypothesis that 
the J /if/ is produced by a hard subprocess where the cc is in a 
colorless state. The dominant hard QCD subprocess for J /if/ 
hadroproduction is thus a 2 — > 2 reaction in contrast to the 
feed-down gg — > x& —> J/^J or CO mechanism such as 
gg — > (cc)g C — > J/i//g [40]. Nuclear shadowing should then 
be implemented along the lines of [47], both for gg and eg 
part, although the c-quark shadowing is poorly known. Thus 
the dedicated study of J /if/ + c in pA collisions could provide a 
unique way to study such shadowing effects as well as heavy- 
quark energy loss. We also note that the yield from eg subpro- 
cesses is expected to have the usual factorizing nuclear depen- 
dence A a(Xl \ where X2 is the light-front momentum fraction of 
the nuclear parton, in contrast to the factorization breaking be- 
havior A a< - XF ^ ~ A 2 ^ 3 observed at high Xp [31, 43], explainable 
by the coalescence of IC pairs turning into CS pairs after in- 
teracting with partons from the target surface [ 44 4 6], 

In conclusion, we have carried out the first NLO analysis 
in the Color-Singlet model of J /if/, if/(2S) and T production at 
RHIC and have shown that the CS yield is in agreement with 
the ^-integrated cross sections measured by the PHENIX 
and STAR collaborations. We have also shown that c-quark- 
gluon fusion is responsible for a significant, and measurable, 
part of the yield, and we call for a dedicated measurement to 
pin down this contribution and assess the importance of the 
charm content of the proton. Such a study may also shed light 
on effects due to color-transfer effects beyond the CSM. We 
predict a significant excess of the lepton yield on the "away" 
side of the J/ifi arising from c-quark jet and argue that the ra- 
pidity dependence of this correlation is strongly sensitive on 
the specific mechanisms for the creation of charm in the pro- 
ton. Finally, we have discussed the implication of our work 
on heavy-ion studies. 
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